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ABSTRACT 

We have investigated the Tully-Fisher relation for Polar Ring Galaxies 
(PRGs), based on near infrared, optical and Hi data available for a sample of 
these peculiar objects. The total K-band luminosity, which mainly comes from 
the central host galaxy, and the measured H l linewidth at 20% of the peak line 
flux density, which traces the potential in the polar plane, place most polar rings 
of the sample far from the Tully-Fisher relation dehned for spiral galaxies, with 
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many PRGs showing larger H l linewidths than expected for the observed K band 
luminosity. This result is conhrmed by a larger sample of objects, based on B- 
band data. This observational evidence may be related to the dark halo shape 
and orientation in these systems, which we study by numerical modeling of PRG 
formation and dynamics: the larger rotation velocities observed in PRGs can be 
explained by a flattened polar halo, aligned with the polar ring. 

Subject headings: galaxies: peculiar - galaxies: kinematics and dynamics 


1. Introduction 

Polar Ring Galaxies (PRGs) are peculiar objects composed of a central spheroidal com¬ 
ponent, the host galaxy, surrounded by an outer ring, made up of gas, stars and dust, which 
orbits nearly perpendicular to the plane of the gas-poor central galaxy (Whitmore et al. 
1990). Previous papers (Arnaboldi et al. 1995; Arnaboldi et al. 1997, lodice et al. 2002a, 
2002b, 2002c) found that even where the morphology of the host galaxy resembles that of 
an early-type system, PRGs show many similarities with late-type galaxies. The PRGs are 
characterised by a large amount of neutral hydrogen (H l), always associated with the polar 
structure (Schechter et al. 1984; van Gorkom et al. 1987; Arnaboldi et al. 1997), and by a 
gas-to-total luminosity ratio in the B-band typical of late-type galaxies. 

The connection between PRGs and spirals is important for two of the best studied 
systems: NGG 660 (van Driel et al. 1995) and NGG 4650A, for which the new surface 
photometry (lodice et al. 2002a; Gallagher et al. 2002), based on near-infrared (NIR) and 
optical Hubble Space Telescope (HST) data, have conhrmed that the polar structure appears 
to be a disk of a very young age. By exploring the properties of the host galaxy and ring 
in the optical and NIR, for a sample of PRGs, lodice et al. (2002a, 2002b, 2002c) found 
that the connection with spirals is tighter. lodice and collaborators found that in almost 
all PRGs i) the host galaxy is bluer (even NIR and optical colors are similar to those of 
late-type galaxies) and younger (age vary between 1 and 5 Gyrs) than normal early-type 
galaxies, and it is characterised by a “compact”, nearly-exponential bulge, with a brighter 
and smaller disk than those found in standard SO galaxies; ii) the polar structure is even 
bluer than the host galaxy, and it has colors similar to those of dwarf irregular and spiral 
galaxies. 

The Tully-Fisher relation (TF) is the most important scaling relation for disks (Tully & 
Fisher, 1977): this is an empirical relationship between the disk rotational velocity (Wot) and 
its absolute luminosity (L), where L oc Wot) approximately. The TF is extensively used to 
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estimate extragalactic distances (e.g. Sakai et al. 2000; Tully & Pierce 2000, and references 
therein). Furthermore, this relation is considered a critical constraint in galaxy formation 
theories (Dalcanton et ah 1997; McGaugh & de Blok 1998; Mo et al. 1998; Steinmetz & 
Navarro 1999, van den Bosch 2000). Several studies have addressed the origin of the TF 
relation (Kauffmann, White & Guiderdoni 1993; Gole et al. 1994; Silk 1997; Avila-Reese, 
Firmani & Hernandez 1998, Heavens & Jiminez 1999; Elizondo et al. 1999; van den Bosch, 
2000), but a dehnitive conclusion is yet to be found. 

In the past few years, several studies have asserted the validity of the TF relation for 
some classes of disk galaxies which show different photometric and kinematical properties 
with respect to ‘classical’, high-surface-brightness spiral galaxies (Matthews, van Driel & 
Gallagher 1998a, 1998b; McGaugh et al. 2000; Ghung et al. 2002). In particular, the studies 
by Matthews et al. (1998b) and McGaugh et al. (2000), while exploring the TF relation 
for extreme late-type spiral galaxies and for faint field galaxies, indicated the validity of 
the “baryonic TF relation”, which accounts for the gas mass in addition to the mass of 
the optical component by making a baryonic correction (Milgrom & Braun 1988). These 
latest developments indicate that the TF relation is probing a very close liaison between the 
dark halo parameters and the total quantity of baryons in galaxies: the dark halo, which is 
responsible for the H l linewidth and the flat rotation curve in the outer regions of a disk, is 
tuned to the total amount of baryons in the luminous component. 

In PRGs, the Hi linewidth (AH) measures the dynamics along the meridian plane, 
which is dominated by the dark matter, while the baryons are nearly equally distributed 
between the host galaxy and the polar ring. We wish to investigate the position of the PRGs 
in the log(AR) — L plane, and study via N-body simulations of 3-D systems whether the 
dark halo shape may influence their position in the log(AI/) — L plane, with respect to the 
TF relation of bright disks. The question of the dark halo shape is important i) to constrain 
dark matter models, through cosmological simulations (Navarro, Frenk & White 1996, 1997; 
Bullock et al. 2001) which predict the distribution of the halo shapes and the universal radial 
dependence of the dark matter distribution; ii) to give hints on the nature of dark matter 
(see Gombes 2002 as a review); and furthermore Hi) the dark halo properties in PRGs can 
give important constraints on the formation scenarios for these peculiar objects, which is 
still an open issue (see lodice et al. 2002a and 2002c and references therein). 

In Section 2 we present the NIR, optical and HI observations available for a sample of 
PRGs and discuss the PRGs location in the log(AR) — L plane; and we derive some insights 
on the dark halo shape in these objects via N-body, 3- and 2-D models in Section 3. In 
Section 4 we evaluate the uncertainties and discuss our Endings with respect to relevant 
dark halo parameters, i.e. total dark mass, central density and flattening , and then derive 



our conclusions in Section 5. 


2. Observations 

2.1. Near-infrared, B band, and Hi observations of polar ring galaxies 

New near-infrared J, H and Kn images are available for a sample of PRGs (lodice et al. 
2002a, 2002b, 2002c), which are selected from the Polar Ring Catalogue (Whitmore et al. 
1990), and for ESO 235-G58, which was classihed as a PRG related object by Buta & Crocker 
(1993). Data were obtained at the 2.3 m telescope of the Mt. Stromlo and Siding Spring 
Observatory, with the CASPIR infrared camera (McGregor 1994). The angular resolution of 
this camera is 0.5 arcsec pixeR^ and it has a held of view of 2.0' x 2.0'. A detailed description 
of the data reduction and analysis is given by lodice et al. (2002a, 2002b, 2002c). The average 
photometric error on K magnitudes is 0.05 mag, and a detailed discussion is given in lodice 
et al. (2002b). This sample contains only 6 PRGs, while more data are available in the 
optical bands. Since we wish to investigate the validity of the TF relation for PRGs, it 
is worthwhile to compare the PRG properties with those of a very large sample of spiral 
galaxies. To derive the TF relation for normal disk galaxies, we will use the very large and 
detailed dataset available in the I-band from Giovanelli et al. (1997). The B band magnitude 
is known for many PRGs (Whitmore et al. 1990; Van Driel et al. 2000a, 2002b; Gallagher 
et al. 2002), and we estimate the B band magnitude for spiral disks in the sample from 
Giovanelli et al. using the observational relation between morphological type index and the 
B — I colors (de Jong 1996). 

In PRGs the bulk of the light comes from the host galaxy, which is a nearly dust free 
component in most systems (see lodice et al. 2002a, 2002b). In most PRGs of our sample, 
the central most luminous part of the host galaxy is not obscured by the dust in the PR 
region which passes in front of it, so there is no signihcant correction for internal extinction. 
For the 3 PRGs (ESO 603-G21, ARP 230, ESO 235-G58) where the dust associated with 
the PR obscures the central regions of the host galaxy, we have corrected the K and B 
total magnitudes for internal extinction, following the equations from Tully et al. (1998). 
We expect little extinction within the polar ring itself, since this component resembles low- 
surface-brightness disk galaxies (see the HST images in lodice et al. 2002a and Gallagher et 
al. 2002), where dust extinction is not high (Matthews and Wood 2001). 

The Hi integrated line prohle data were obtained from several published hydrogen 
observations of PRGs, carried out by e.g. Richter et al. (1994), van Gorkom et al. (1987), 
van Driel and collaborators (2000, 2002a, 2002b), with several radio telescopes. 
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The absolute Kn magnitude {MKn), B band absolute magnitude, and measured linewidth 
at 20% of the peak line flux density (AV 20 ), foi' each PRG in our sample are listed in Table 1. 
For all distance-dependent quantities we have assumed Hq = 75 km s“^ Mpc“^. 


2.2. PRGs and the TF relation for spiral galaxies 

We wish to compare the position of the PRGs in the log(AR) — L plane with the 
TF relation observed for bright spiral galaxies. In principle, we do not know what expect, 
because the two quantities, rotational velocity and absolute luminosity, are related to two 
different components in the PRG, the polar structure and the host galaxy respectively, which 
he in two decoupled planes. In fact the total K light comes mostly from the central galaxy 
(lodice et al. 2002a, 2002b, 2002c), whereas the observed Hi linewidth traces the dynamics 
of the polar structure along the host galaxy’s meridian plane. 

In Figure 1 we show the K-band TF relation for a sample of spiral galaxies studied by 
Verheijen (1997, 2001). The values of Mk and log(AV 2 o) for PRGs are also shown on this 
plot. Our and Verheijen’s data sets have very similar photometric properties and limiting 
magnitudes. 

We see that hve PRGs he near to the high-velocity boundary of the TF relation, or show 
larger velocities (for a given luminosity) than disk galaxies. Only the PRG AM2020-504 
shows a lower rotation velocity for its K-band absolute magnitude. 


In the B-band, the tendency of PRGs to have larger velocities with respect to the TF is 
conhrmed for a larger sample of PRGs compared with data for 787 disk galaxies (Giovanelli 
et al. 1997). In Figure 2, we see that two PRGs lie at lower velocities than those predicted 
by the TF for spiral galaxies, two objects lie on the TF relation, and twelve PRGs either 
lie on the high velocity boundary of the TF relation or show much larger velocities. Van 
Driel et al. (2002b) also hnd that most kinematically conhrmed PRGs show larger Hi prohle 
widths than bright spiral galaxies, at a given luminosity. 


3. The PRG positions in the log(AV) — L plane and their implications for the 

dark halo 

Observing higher or lower velocities with respect to the linear TF of disk galaxies is 
relevant for the discussion on the dark halo shape. Via analytical models and simple as¬ 
sumptions about the mass distribution, either luminous or dark, we can estimate where the 
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PRGs ought to lie in the log(AR) — M^/b plane (Fig.l and Fig.2), with respect to the TF 
relation for disk galaxies. If there were no dark matter, and the gravitational potential were 
oblate in the same sense as the flattened host galaxy, the polar ring would acquire an eccen¬ 
tric shape. When the polar ring and the host galaxy are both seen edge-on, which is close to 
being the case for most of our PRGs, the net effect will be that the line-of-sight (LOS) polar 
ring velocities are reduced, see Fig. 3. In the logarithmic, scale-free, potential case, a simple 
formula gives the expected velocity ratio between the major and the minor axis components 
as: 

= 1- — = ep ~ 2e<i, (1) 

'^minor 

from Gerhard & Vietri (1986), where is the flattening (1-axis ratio) of the density distribu¬ 
tion and is the potential flattening. Therefore we would expect PRGs to have on average 
lower velocities with respect to what would be measured in the equatorial plane: see Fig.4. 
This implies that when the polar structure is eccentric, the observed LOS velocities in AV 20 
are the smallest, i.e. those from the particles in the polar regions, as shown in Figure 4. 
Thus the observed AV 20 depends on both the mean velocity along the ring, and the ring 
eccentricity. On the contrary, Fig.l and Fig.2 show that the majority of PRGs have larger 
velocities than expected in the log(AR) — Mk/b plane. Therefore we need to investigate 
how these velocities can be produced, and how they may depend on the intrinsic properties 
of the dark galaxy halo. 

We hrst compute a series of N-body models of the formation of polar rings in galaxy 
mergers or by tidal accretion of gas from another galaxy (Bournaud & Gombes 2002). Both 
scenarios are known to form polar structures (see Reshetnikov & Sotnikova 1997; Bekki 1997, 
1998). The numerical code computes the gravitational potential via an FFT method (James 
1977), and includes sticky particles gas dynamics (Schwarz 1981), and star formation that is 
assumed to obey a generalised Schmidt law (Schmidt 1959). Both scenarios for formation of 
polar rings give similar results for the kinematics of the ring after its formation. The details 
of the N-body models and numerical codes are described in Bournaud & Gombes (2002). 
Fig.5 shows an example of the structures produced. 

These 3D N-body models do not allow us to choose the dark halo shape after the polar 
ring is formed. Moreover, we cannot separate the influence of the host galaxy potential, the 
dark halo flattening, and the ring self-gravity, which are all important in order to understand 
what our observations may imply for the PRG structure. We therefore developed a simpler 
model, where every parameter can be varied. In this model, the dark halo and the host 
galaxy are two rigid potentials. The eccentricity of the polar structure is given by the shape 
of the particle orbits in this rigid potential. Particles are launched around the polar plane 
and are followed for 4 Gyrs, before their kinematics is analysed. The polar ring self-gravity 
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may also be included via a 2-D FFT algorithm that computes the gravitational potential in 
the polar plane. 

Our N-body 3D models of polar rings (see Figure 5) show different characteristics de¬ 
pending on the shape of their dark halo. We see that when the halo is oblate and flattened 
towards the host galaxy, the observed velocities are then smaller, and the PRGs lie on the 
left side of the TF relation of high-surface-brightness disks. When the oblate halo is flattened 
towards the polar ring, the observed velocities are larger, shifting the PRGs to the right side 
of the diagram. As most of the observed PRGs have larger velocities compared to the TF 
relation for spiral galaxies, it seems to suggest that their dark halo is flattened towards the 
polar plane. 

The second model is initially computed for massless rings. We choose a lOkpc-scaled 
host galaxy, and a polar ring whose mean radius is 15 kpc. The dark halo is computed 
to produce a flat rotation curve at large radii, its mass is about 0.8 times the host mass 
inside 10 kpc radius, and 1.4 times the host galaxy mass inside 15 kpc. A halo flattened 
towards the equatorial plane reduces the measured spread in radial velocity in the polar 
structure: the more flattened the potential, the more eccentric the ring is and the smaller 
the observed velocities are. When the halo is spherical, the potential is flattened by the host 
galaxy luminous mass, therefore the ring is still eccentric, and the PRG velocities then he 
on the low-velocity boundary of the TF relation, far from the observational result. This is 
what was found by Gombes and Arnaboldi (1996) in their detailed model for the polar ring 
galaxy NGG4650A: a spherical halo cannot reproduce the large ring velocities, because the 
ring is eccentric, unless when the halo is very massive (but this case is easily ruled out, see 
Sect. 4.4). 

We now wish to account for the ring self-gravity: the ring mass for all objects in our 
sample is believed to be a significant fraction of the stellar mass of the host galaxy or of 
the same order, as it is in the case for NGG 4650A (lodice et al. 2002a). When the halo 
is spherical or flattened towards the host galaxy equatorial plane, its self-gravity makes the 
polar ring more circular (see Figure 6). In Figure 1 we have shown the effects for polar rings 
whose mass is equal to the mass of the host galaxy. When the polar ring is lighter, these 
effects are smaller; for instance when the polar mass is half the mass of the luminous HG, we 
hnd results that are comparable to those obtained with the first model. When self-gravity 
is considered, the simulated PRGs he near the TF relation for bright disk systems of Fig.l. 
And yet, when the flattened dark halo is aligned with the host galaxy equatorial plane, the 
ring’s self-gravity cannot account for the large velocities observed in most polar rings, even 
when the ring is as massive as the host galaxy. But, when dark halos are flattened towards 
the polar structures, the ring self-gravity is less important, because the flattened potential 



already makes rings more circular: the positions of the corresponding points does not move 
much from the the massless ring case, when the ring mass is increased, see Figure 1. Indeed 
such systems seem able to explain the large velocities observed in the polar structures, and 
conhrms the conclusions drawn from our 3-D accretion and merging models. 

Let us now consider the case of AM 2020-504 and UGC 4261, which show lower velocities 
than bright disks, for their observed K and B band absolute magnitudes. Both galaxies have 
narrow rings with very small masses, therefore they are well described by our massless ring 
model. This model easily predicts the observed velocities in the polar ring, when the potential 
is signihcantly flattened towards the host galaxy equatorial plane. This is conhrmed in the 
case of AM 2020-504 by the photometric study of the host galaxy, which is a flattened 
elliptical (E4, see Arnaboldi et al. 1993). Moreover, if those systems are embedded in 
massive dark halos that contribute much to the potential near the polar ring, these halos 
must be flattened too. 


4. Discussion 

In the log(AD) — M^/b planes shown in Fig. 1 and Fig. 2, most PRGs have larger Hi 
rotation velocities than standard spiral galaxies, at a given K or B-band luminosity of the 
stellar component. Our N-body simulations have suggested that a likely explanation for this 
effect is a flat dark halo, whose main plane is aligned with the host galaxy meridian plane, 
and prevents the polar ring to become eccentric. The question arises if other effects, i.e. 
non-homogeneities in the TF relations for spirals, caused by bar and/or non edge on disks, 
or larger M/L ratios, can produce similar results and therefore be alternative explanations 
for the high velocities observed in PRGs. 


4.1. Possible biases in the spiral sample 

PRGs are peculiar objects that are easily identihed when both their central component 
and polar ring are observed nearly edge-on. When dealing with the TF relation for spiral 
galaxies, we have assumed that it was valid for axisymmetric edge-on disks. Yet, the large 
sample from Giovanelli et al. (1997) contains barred disks and disks at lower inclination 
with respect to the LOS. Bars may cause a bias in the observed velocities, and calculating 
the real value of AV for a disk at a lower inclination cannot be done exactly. We have 
therefore selected from the Giovanelli et al. sample those disk galaxies that are not barred 
and are seen nearly edge-on [i > 80); we built a new TF relation from this new sample of 
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179 disks, which is now free from those biases that may cause lower Hi velocities. The linear 
interpolation of this new TF relation is shown in Figure 2 : it is still very close to the one of 
the original sample and indicates that barred disks and disks with lower inclination do not 
modify the mean TF relation (although the scatter is smaller, for bars increase the scatter 
of the TF relation). So this cannot explain why PRGs have larger observed Hi velocities 
than bright disks. 


4.2. The baryonic Tully-Fisher relation 

We have already mentioned (Sec. 1) that PRGs are gas-rich, thus we wish to investigate 
whether once we account for the gas mass, this might move them back on the TF relation, 
i.e. if they are close to the baryonic TF relation (Matthews et ah 1998b; McGaugh et al. 
2000). A mass-to-velocity TF relation for PRGs is shown in Figure 7: PRGs are now closer 
to the TF relation for disk galaxies, but still do not lie on this relation: i.e. the higher 
velocities in the PR cannot be explain by additional mass in the Hi associated within the 
PR itself. 


4.3. Host galcixy mass-over-light ratios 

Another possible explanation for higher H l rotation velocities is that PRGs have larger 
mass-to-light ratios than spiral disks of similar K&B-band luminosity, i.e. for the same 
observed AV 20 , PRGs may be less luminous on average in K/B than standard spiral disks. 
Therefore we have computed the linear interpolation of the TF relation for the Sa galaxies 
alone, of the Giovanelli et al. (1997) sample. The average B-band magnitude of this sample 
corresponding to log(AR) = 2.5 is -19.35, while it is brighter, -19.16, for the whole Giovanelli 
sample. Still, the difference of 0.19 mags is not enough to place most PRGs on the TF 
relation, which would require a value 5 times larger. 


4 . 3 . 1 . Comparison between PRG host galaxies and SOs 

The similarities between PRGs and spiral galaxies, from NIR and optical photometry, 
prompted us to compare PRGs with the TF for early-type spiral galaxies rather than with 
lenticulars. But we have also compared the position of PRGs in the log(AR) — Mb plane 
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with the TF relation for SO galaxies^ (Neistein et ah 1999, Mathieu et ah 2002), and we 
hnd that many PRGs lie near the TF relation for SO galaxies. The offset between the TF 
relation for SO galaxies and the TF for bright spiral galaxies is caused by the SOs being less 
luminous than bright spirals, at a given velocity. As suggested by Mathieu et al. (2002), this 
is consistent with the formation scenario in which SO galaxies were formed from ordinary 
late-type spirals which were stripped of their star-forming medium. 

Can this be a valid explanation for the PRGs position in the log(AR) — Mk/b plane? If 
the PRG positions in this plane were caused by larger M/L ratios for the luminous compo¬ 
nent, as is the case for lenticular galaxies, it must be so for the host galaxy too, which is the 
component most similar to an SO galaxy. We have computed log(AR) for the host galaxies 
of the best studied systems, NGC 4650A, NGC 660, NGC 2685, UGC 7576 and A0136-0801, 
(Sackett et al. 1994; van Driel et al. 1995; Simien & Prugniel 1997; Whitmore et al. 1990; 
Schweizer et al. 1983), using optical absorption-line rotation curves along the host galaxy 
equatorial plane; these are shown in Fig. 8. The host galaxies in these hve systems fall on 
the TF relation of bright spiral disks, which indicates that i) the mass-to-light ratio of this 
component is different from those of standard SOs, as also found in recent studies of PRG 
NIR and optical colors (see lodice et al. 2002a, 2002b, 2002c) and ii) a luminous-to-dark 
matter content similar to those of standard bright disks. 


4.4. Peculiar properties of the dark halo 

Gerhard et al. (2001) showed that elliptical galaxies follow a TF relation in the log(AI/) — 
M*, where M* is the total mass in the luminous component, which is shallower than the re¬ 
lation for spiral galaxies, even when the maximal M/Lb is adopted to compute the total 
stellar masses. This led Gerhard and collaborators (2001) to infer that elliptical galaxies 
have slightly lower baryonic mass than spiral galaxies of the same circular velocities, and 
that their dark halos are denser than halos of spiral galaxy with the same Lb- How much 
more massive must the dark halo be, to account for the velocities observed in polar rings? 
In our numerical models, we had assumed a dark-to-visible mass ratio of about 1 inside 
the optical radius^, as it is the case in spiral galaxies, and only flattened polar halos could 
reproduce the observed Hi velocities. The observed value of AV 20 depends largely on the 
dark halo shape and the ring eccentricity, while it varies only as the square root of the total 


^The TF relation for SO galaxies was obtained by comparing the asyntotic rotation speeds of stars with 
the absolute magnitude (see Mathieu et al. 2002). 

^For PRGs, the optical radius correspond to the ring average extension. 
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mass (and depends even less on the dark mass). Thus, a large amount of dark matter is 
needed when the halo is not polar: if one wishes to reproduce the PRG position on the TF 
diagram with a spherical halo, then the dark mass needed is twice than the visible mass. 
With an E5 equatorial halo, the dark-to-visible mass ratio raises to about 3.5. 

An increased total dark halo mass alone can explain the observed position of PRGs on the 
log(AR) — Mb plane, but the required large dark masses are unrealistic: the fact that the 
PRGs host galaxies lie within the TF for spiral galaxies (see Fig. 8) indicates that the dark 
matter content in PRGs is similar to that of bright spiral galaxies. 

Also in the spherical halo case, the host galaxy would be offset too from the TF relation 
of bright disks, because the dark mass. In the case of NGG 4650A, Gombes & Arnaboldi 
(1996) found that no dark matter was required to reproduce the host kinematics, and that 
the dark-to-visible mass ratio could not be as large as required by a spherical halo. If we 
want to account for observed velocities in PRGs with a spherical halo, we have to believe 
that the halo density drops more slowly than 1/r. 

A prolate dark halo associated with the host galaxy, i.e. with the major axis aligned 
with the HG major axis, could also account for the large velocities observed in rings, provided 
it is massive enough. To reproduce the large velocities in the PR, such a halo must be even 
more massive than for the oblate case. For an E5 prolate halo and a ring mass that amounts 
to 25 % of the visible mass, we hnd that the dark-to-visible mass ratio inside the optical 
radius would be close to 4. Such a massive halo would cause a large offset of the host galaxy 
from the TF relation of bright disks, which is not observed. On the other hand, a prolate 
halo with a polar major axis can account for the velocities in PRGs and this is one of the 
possibility of a halo flattened towards the polar ring. 

Given the constraint from the host galaxy observed kinematics, the dark halo shape 
alone can explain the high rotation velocities observed in polar structures, provided that 
the short axis of the halo lies perpendicular to the polar ring, in the equatorial plane of the 
central galaxy. In the case of NGG 4650A, Gombes & Arnaboldi (1996) had already argued 
that a dark halo flattened towards the polar ring required less dark mass to account for 
the large velocities observed in the polar ring, while the visible mass alone was enough to 
account for the velocities in the host galaxy. 


5. Conclusion 


In this paper we have explored the position of PRGs in the log(AR) — M^/b plane 
with respect to the TF relation for spiral galaxies. The majority of PRGs show larger 
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rotation velocities at a given luminosity than spiral galaxies. A comparison with the N-body 
simulations has led us to conclude that these large velocities indicate that the dark halo is 
most likely flattened towards the polar ring plane. This explains the position of PRGs with 
respect to the TF relation for normal spiral disks, the position of the host galaxy within 
the TF of bright disks, and it may become an important constraint on the nature of dark 
matter. 


The 3-D simulations in this work were computed on the Fujitsu NEC-SX5 of the CNRS 
computing center, at IDRIS. This research made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and Space Administration. M.A. 
wishes to thank O. Gerhard for many useful discussions on the TF relation for elliptical 
galaxies and for his comments on an earlier version of this paper. 
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log(AV2o) 

Fig. 1.— Absolute magnitude in the K band vs. the measured Hi linewidth at 20% of the 
peak line flux density (AV 2 o); foi’ PRGs of the selected sample, compared with a sample 
of spiral galaxies from Verheijen (2001), and with the results from 3-D simulations and 2- 
D models (for massless rings and rings that are as massive as the host galaxy). We have 
assumed Hq = 75 km s“^ Mpc“^, as made also by Verheijen. The long-dashed line is a linear 
interpolation of the TF relation for spiral galaxies, and the short-dashed lines show the width 
at 15% of the peak of the statistical distribution of spiral galaxies. For both models, the 
flattening of the halo is indicated next to each circle (massless ring) or triangle (very massive 
ring): a positive x number indicates that it is an Ex halo with an equatorial flattening, while 
—X corresponds to an Ex halo flattened toward the polar plane. The results from the 3-D 
models shown in this plot are those computed for the accretion scenario (Reshetnikov & 
Sotnikova, 1997) ; our values for AV 20 vs. Mi^- are very similar when one considers the 
merging scenario (Bekki, 1998). 
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Fig. 2.— Absolute magnitude in the B band vs. the linewidth at 20% of the peak line 
flux density (AV 20 ), for 15 PRGs. Data for disk galaxies (dots) are from Giovanelli et al. 
(1997). Absolute magnitudes have been normalised to the same value of Hq for PRGs and 
disk galaxies (75 km s“^ Mpc“^). A linear interpolation of the TF relation is shown for these 
disk galaxies (solid line), 81% of which lie inside the dashed lines that are computed at 25% 
of the peak of the statistical distribution of spiral galaxies. The long-short dashed line is 
obtained for unbarred disks, seen nearly edge-on (i > 80, see Section 4). 







18 


line of sight 


observed velocities 


host galaxy 


cr ^ 


^polar ring 


Fig. 3.— Line of sight for most PRGs: both components are seen nearly edge-on, thus the 
measured velocities in the polar structure are the smallest ones when the ring is eccentric, 
i.e. when the potential is flattened along the equatorial plane. This line of sight corresponds 
to the PRGs of our sample, and it is adopted for the analysis of the numerical models. 
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Fig. 4.— Simulated Hi profiles for circular and eccentric polar rings seen edge-on, as shown 
in Figure 3. Solid line: Hi prohle computed for a circular ring, extending from 10 to 15 
kpc. Dashed line: Hi profile computed for an eccentric ring, with the same radial extension 
and radius, and ellipticity 0.35 (similar to the ring from 3-D simulations shown in Figure 5). 
For gas, we assume a velocity dispersion of 10 km/s. The values of log(AV 2 o) are 2.51 for 
the circular ring and 2.38 for the eccentric ring. This difference is in agreement with the set 
of values for AV 20 , at a given total K or I band luminosity, obtained from PRG numerical 
models (see Figure 1). 
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Fig. 5.— Example of polar ring formed in the 3-D N-body model, according to an accre¬ 
tion scenario. In most observed PRGs both the central disk and the polar ring are seen 
nearly edge-on (right panel). This polar ring is eccentric, since the gravitational potential is 
flattened by the host disk and an E4 dark halo. 
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Fig. 6.— Polar rings from our 2-D model, formed during an accretion event, obtained with 
a E3 halo flattened in the equatorial plane {x — y). The ring lies in the z — x plane. The 
massless ring (on the left) is eccentric, which reduces the observed velocities. Though the 
very massive ring (on the right, the ring is as massive as the host galaxy) is rounder, its self¬ 
gravity is not sufficient to reproduce the observed velocities, unless the halo is also flattened 
towards the ring. The measured linewidth AV 20 is 30% larger for the massive polar ring 
than for the massless ring. 



Fig. 7.— Plot of the stellar, mass in solar units, (smaller symbols for PRGs, long-short 
dashed line for spirals) and the stellar plus gaseous mass (solid line and dots for spirals, 
large symbols for PRGs) against the observed rotation velocities. Due to the similarities 
between host galaxies in PRGs and early-type spiral galaxies (Arnaboldi et ah 1995; lodice 
et al. 2002a, 2002b), we have assumed M/L = 2 for this component, as for spirals (de Jong 
1996). In this plot, data for disk galaxies are from Mathewson & Ford (1996). 
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Fig. 8.— TF relation for SO galaxies (data from Neistein et al. 1999 and Mathieu et al. 
2002) and PRGs. The TF for spiral galaxies is plotted as a dashed line. Large crosses show 
the position of the central component in hve PRGs. For NGG 4650A the circular velocity 
was derived from the rotation velocities and velocity dispersion of the stellar component 
(Gombes & Arnaboldi 1996), by computing the corresponding gas linewidth and assuming 
a 10 km/s dispersion (which is the standard value for disks). The total B magnitude for 
NGG 4650A is from lodice et al. 2002a. For NGG 660, velocities are derived from Ha data 
(see van Driel et al. 1995). For NGG2586 and A0136-0801 the B magnitudes relative to the 
HG are given by Reshetnikov et al. 1994. The arrow shows the offset in log(AV 2 o) between 
the velocities in the host galaxy equatorial plane and in the polar ring for each system. 





Table 1. 


Absolute Kn and B band magnitudes, heliocentric velocity, and Hi linewidth at 
20% of the peak line flux density for a sample of polar ring galaxies. 


Object name 

PRC 

V 


Ms 

AV20 

Ref. 



krn.s”^ 

mag 

mag 

km.s 



-19.23 365 vD 


A0136-0801 

A-01 

5528 

ES0415-G26 

A-02 

4604 

NGC2685 

A-03 

871 

UGC7576 

A-04 

7022 

NGC4650A 

A-05 

2910 

UGC9796 

A-06 

5420 

ARP230 

B-01 

1710 

IG1689 

B-03 

4557 

ESO503-G17 

B-12 

10481 

UGC9562 

B-17 

1242 

AM2020-504 

B-19 

4963 

ESO603-G21 

B-21 

3124 

NGC660 

C-13 

850 

UGC4261 

C-24 

6408 

UGC4385 

C-27 

1969 

ES0235-G58 


4262 


-22.69 

-19.24 

360 

vD 


-18.20 

335 

vD 


-19.02 

465 

vD 

-21.37 

-18.67 

235 

vD 


-18.25 

355 

vD 

-20.95 

-18.40 

245 

vD 


-20.40 

300 

vG 


-19.19 

415 

vD 


-16.68 

200 

vD 

-22.63 

-19.60 

182 

vD,A 

-21.74 

-18.60 

285 

vD 


-18.25 

325 

vD 


-19.87 

240 

vD 


-17.58 

195 

vD 

-22.13 

-19.21 

292 

vD 


^References for values of AV 20 • telescope codes (V:VLA, N: Nangay, G: 
Green Bank, P: Parkes); vG: van Gorkom et al. (1987), vD: van Driel et 
al. (2000, 2002a, 2002b and references therein) and Richter et al. (1994), A: 
Arnaboldi et al. (1993). 

^The B band magnitude for NGG 4650A is from Gallagher et al. (2002). 



